#!/usr/bin/gnuplot
set encoding iso_8859_1
set term pdf enhanced color solid size 4in, 4in

#============================
#F(Q_LC), F(Q_NC), F(Q_LN)
#============================
#set xlabel "Time Steps"
#set x2label "Free Energy (K_bT)"
#set x2tics
#set ylabel "Q_{Aspe}"
#set output "FQAspe.pdf"
#p 'QPLU.smth' u 1:2 w l ti "Q_{Aspe}",#  'QPLU.smth.wam2' u 2:1 w l axis x2y1 lw 4 lt 7 ti "F({Q_{Aspe}})"
#
#set ylabel "Q_{Ospe}"
#set output "FQOspe.pdf"
#p 'QPLU.smth' u 1:3 w l ti "Q_{Ospe}",#  'QPLU.smth.wam3' u 2:1 w l axis x2y1 lw 4 lt 7 ti "F({Q_{Ospe}})"
#
#set ylabel "Q_{Hspe}"
#set output "FQHspe.pdf"
#p 'QPLU.smth' u 1:4 w l ti "Q_{Hspe}",#  'QPLU.smth.wam4' u 2:1 w l axis x2y1 lw 4 lt 7 ti "F({Q_{Hspe}})"
#
#set ylabel "Q_{C0L2H}"
#set output "FQC0L2H.pdf"
#p 'QPLU.smth' u 1:5 w l ti "Q_{C0L2H}",#  'QPLU.smth.wam5' u 2:1 w l axis x2y1 lw 4 lt 7 ti "F({Q_{C0L2H}})"
#
#============================
#F(Q_LC,Q_NC), F(Q_LC,Q_ATP), F(Q_NC,Q_AMP)
#============================
set pm3d at bs
set palette color
#set palette model CMY
set palette defined ( 0 "white", 1 "green", 2 "cyan", 3 "blue")
set view map
#set cblabel 'Energy Bar' textcolor lt 3

set title "CA Contact Map"
set xlabel "Protein Residue Index"
set ylabel "RNA Residue Index"
set xrange [0:72]
set yrange [71:114]
set output "contactmap.pdf"
splot "nonnative.dat" using 4:5:3 w p pt 6 ps 0.8 palette noti,\
      "native.dat" using 4:5:6 w p pt 5 ps 0.3  palette noti

